Probing the pairing symmetry of the iron pnictides with electronic Raman scattering 
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An important issue in the study of the iron-arsenic based superconductors is the symmetry of 
the superconducting gap, a problem complicated by multiple gaps on different Fermi surface sheets. 
Electronic Raman scattering is a flexible bulk probe which allows one in principle to determine gap 
magnitudes and test for gap nodes in different regions of the Brillouin zone by employing different 
photon polarization states. Here we calculate the clean Raman intensity for A\ g , B\ g and B29 
polarizations, and discuss the peak structures and low-energy power laws which might be expected 
for several popular models of the superconducting gap in these systems. 
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I. INTRODUCTION 

Since their discovery 1 , there has been a considerable 
effort to understand the origin and nature of supercon- 
ductivity in iron-pnictide materials (for early reviews see 
Refs. 2-4). Initial information on the structure of the 
gap is often provided by power laws in the temperature 
dependence of thermodynamic and transport properties, 
which are related to the topology of the superconduct- 
ing gap in its nodal regions. Nuclear magnetic resonance 
(NMR) studies 5-8 showed a low temperature T 3 spin lat- 
tice relaxation rate typical of a gap with nodes. However, 
penetration depth measurements 9-15 have been fit both 
to exponential activated T-dependence, indicative of a 
fully gapped state, and low-T power laws. ARPES mea- 
surements on single crystals of 122-type materials 16-21 
measured the spectral gap reporting isotropic or nearly 
isotropic gaps on all Fermi surface sheets. It is possible 
that these differences reflect genuinely different ground 
states in different materials. However, the complex in- 
terplay of multiband effects, unconventional pairing, and 
disorder leaves open the possibility that a single ground 
state symmetry exists in all the Fe-pnictides, and that 
differences in measured properties can be understood 
by accounting for electronic structure differences and 
disorder 22 ' 23 . It is likely that a consensus will be reached 
only after careful measurements using various probes on 
the same material, and systematic disorder studies. 

We argue here that measurement of the electronic Ra- 
man scattering in the superconducting state can provide 
important information on the structure of the bulk su- 
perconducting order parameter through its sensitivity to 
both symmetry and gap scales. Here we discuss the Ra- 
man scattering for some simple models of the Fe-pnictide 
superconductors. In general the energy of the peaks in 
the Raman intensity are directly related to the magni- 
tude of the gaps on the various Fermi sheets. However, 
whether a given Raman polarization weights a given gap 
strongly or weakly depends on the polarization state of 
the measurement via the Raman vertex 7k. This is par- 



ticularly important for superconductors where the gap is 
strongly momentum dependent, and was exploited suc- 
cessfully in cuprate materials to help determine the d- 
wave symmetry in those systems 24-27 . In addition, the 
presence of nodes and the dimensionality of the nodal 
manifold may be determined by comparison with low en- 
ergy power laws in the Raman intensity in different po- 
larization states. 

In the Fe-pnictide materials, density functional theory 
calculations 28 ' 29 for the paramagnetic state have gener- 
ally found a multisheeted Fermi surface consisting of con- 
centric, nearly circular hole pockets around the T point 
(here referred to as "a sheets" ) and nearly circular elec- 
tron pockets around the M points, or X point in the 
unfolded, 1-Fe zone (referred to as "/3 sheets"). These 
sheets are plotted in Fig. 1. Spin fluctuation mod- 
els of pairing 30-33 have usually found that the leading 
pairing instability had s-wave symmetry, and noted that 
the next leading channel had d x 2_ y 2 symmetry. Wang 
et al. 34 , who studied the pairing problem using a func- 
tional renormalization group approach within a 5-orbital 
framework, also found that the leading pairing instability 
occured in the s-wave channel, and that the next leading 
channel had d x 2_ y 2 symmetry. These calculations differ 
according to whether or not actual nodes are found in the 
leading pairing states, but they all find large anisotropies 
over the Fermi surface sheets, not anticipated in the orig- 
inal predictions of extended-s-type states with isotropic 
gaps which changed sign between the a and (3 sheets 35,36 . 

While ARPES experiments 16-21 have reported weakly 
anisotropic gaps, this may be due to momentum reso- 
lution issues and the difficulty of preparing good sur- 
faces in these systems at present. In particular, the sam- 
ples used in these experiments may be sufficiently dirty 
at the surface that considerable momentum averaging, 
with concomitant gap averaging, could be taking place 23 . 
Angle-dependent specific heat measurements in a mag- 
netic field 37 , which probe the bulk, may be feasible in 
the future but are difficult at present due to the require- 
ment of large clean crystals. Raman scattering with use 
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that due to the form for the Raman vertices, which are 
allowed when there arc multiple Fermi surfaces, there can 
also be backflow effects on symmetry channels other than 
A\ g in the Fe-pnictides. 

We begin in Section II by considering model one-band 
clean systems with gaps inspired by proposals for the 
Fe-pnictides to illustrate what intuition we can gain re- 
garding the Raman response for various polarizations. 
In Section III, we consider gaps on all four fermi surfaces 
predicted by density functional theory. We present our 
conclusions in Section IV. 



II. ELECTRONIC RAMAN SCATTERING IN 
CLEAN SYSTEM 

A. General Theory 



FIG. 1: Fermi surface of model Fe-pnictide represented in un- 
folded (1-Fe) Brillouin zone. The crystalline a, b axes are in- 
dicated in blue, and the polarization geometries for incoming 
and scattering photon polarizations are denoted for Bi g and 
E>2 g geometries. Note that our symmetry notation is rotated 
by 45 degrees with respect to the lattice symmetries. 



of different polarizations may therefore be the best cur- 
rent method of acquiring momentum-dependent informa- 
tion on the structure of the bulk superconducting order 
parameter. We discuss several cases below which should 
allow extraction of the crude momentum dependence and 
possible nodal structure of the order parameter over the 
Fermi surface. Some of these cases have been considered 
in an earlier paper by Chubukov et al. 38 , who however 
focussed solely on the A\ g polarization, and examined 
vertex corrections due to short-range interactions. How- 
ever, they neglected "charge backflow" effect of the long- 
range part of the Coulomb interaction which is required 
to ensure charge conservation. Here we show that for the 
A\ g case, significant changes are to be expected due to 
the coupling of the A\ g Raman charge fluctuations due 
to the backflow effects. 

While the expressions for the B\ g and B^ g channels 
are generally well-characterized by the bare bubble cal- 
culation for the cuprates, it is well-known that the A\ g 
contribution is significantly more complicated due to the 
issue of charge backflow as noted above, and to the num- 
ber of different excitons which may be pulled down from a 
gap edge condensate. In systems with several condensate 
pairing instabilities lying in close proximity in parame- 
ter space, strong excitonic peaks may occur. While the 
possibility of such excitonic modes is of interest for the 
pnictides, in this paper we focus on the generic quasi- 
particle features for Raman scattering, and will consider 
excitonic peaks in a future publication. Charge back- 
flow and Coulomb screening is explicitly included in all 
A\ g calculations, but pair interaction corrections will be 
neglected. We will however also consider the possibility 



Raman scattering is the inelastic scattering of polar- 
ized light from a material. For a review see Ref. 24. 
The cross section of the scattered light is proportional 
to the imaginary part of the channel-dependent Raman 
susceptibility 
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Here we will take a simple frequency independent form 
for the Raman vertices (non-resonant scattering) and 
write the effective Raman charge fluctuations in the 7- 
channel as 
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(k) defines the 

momentum- and polarization-dependent Raman vertices, 
which may include intra- and interband transitions. Gen- 
erally, the vertex is determined by both density and cur- 
rent matrix elements between the conduction band and 
the excited states, and has not been calculated for even 
simple metals like Al. However, the polarization geome- 
tries of the incoming and outgoing photons impose an 
overall symmetry due to the way in which excitations 
are created in directions determined by the electric field 
oscillations, and classifications of the anisotropy of the 
Raman vertices can be employed. 

In this paper, in order to focus on general features for 
Raman scattering in the pnictides, we will neglect band 
structure features and treat all Fermi surface sheets as 
circles. This allows for a simple symmetry classifications 
for the Raman vertices as has been done in the cuprates. 
Expanding the polarization-dependent vertices in Fermi 
surface harmonics for cylindrical Fermi surfaces, 



7n(<?) Al9 = a n + b n cos(46<) 



ln(0) 



c n cos(2#) 
d n sin(20), 



(3) 



3 



with anglc-indcpcndcnt band prefactors a ni b ni c ni d n set- 
ting the overall strength of the Raman amplitudes for 
band n. Since an isotropic density fluctuation vanishes 
for q — ► 0, it can be shown via Eq. 4 that the A\ g con- 
tribution a n to the Raman vertex is cancelled by charge 
backflow and does not contribute to the scattering cross- 
section. This leaves the cos(4#) as the first non-vanishing 
contribution in the expansion. 

As shown in Fig. 1, we note that our choice of coordi- 
nates rests on the 1 Fe unit cell, which is rotated by 45 
degrees with respect to the 2 Fe unit cell. Thus our sym- 
metry notation is electronic and not associated with the 
lattice principal directions, and therefore our classifica- 
tions are 45 degrees rotated with respect to conventional 
lattice classifications. Using lattice coordinates, what we 
call Big would be B 2g , and vice- versa. While this might 
create some confusion, it is convenient to understand the 
interplay of the angular dependence of the vertices and 
the energy gaps in the rotated 1 Fe unit cell frame, as 
shown in Fig. 1. This should be kept in mind however 
when one discusses, e.g., electronic excitations together 
with lattice excitations. Then our symmetry labels B\ g 
and Big would have to be interchanged. 

Other forms for the Raman vertices are allowable, 
with the only requirement being that they must obey 
the transformation rules according to the relevant point 
group symmetries of the crystal. We note that for multi- 
sheeted Fermi surfaces shown in Fig. 1, there are different 
possible forms for the vertices other than Eq. 9. While 
the Raman vertices of the a sheets must transform ac- 
cording to Eq. 9, for the (3 sheets other forms for the 
vertices could be admissable. This includes, for the B lg 
vertex for example, a vertex which is momentum inde- 
pendent on each (3 sheet but of opposite sign. For the 
Big vertex, a Raman vertex which is p-wave like on each 
(3 sheet (with a change of sign) would also be admissable. 
We will explore these possibilities in Sec. III. 

For n bands crossing the Fermi level, the intraband 
Raman response in the absence of Coulomb screening 
and charge backflow is given by 
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is the Tsuneto function for the n th band, having band 
dispersion e„(k), energy gap A„(k), and quasiparticlc 
energy E^(k) = e 2 (k) + A 2 (k). Taking the imaginary 
part of Eq. (5) we then obtain for the Raman response 
at T = 
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Since Raman scattering probes charge fluctuations in 
the long- wavelength limit, the role of the long-range 
Coulomb interaction is important. Isotropic charge fluc- 
tuations are coupled across all unit cells, and Raman 
scattering at low energies must vanish due to particle 
number conservation, leaving only an inelastic light scat- 
tering peak at the plasmon energy. Screening by the long- 
range Coulomb interaction can be taken into account by 
including couplings of the Raman charge density p to the 
isotropic density p fluctuations and is given by 
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Eqs. (4-9) constitute closed form expressions for the 
intraband, non-resonant contribution to the Raman re- 
sponse. 

It is clear that the Raman response is in general not 
simply additive with respect to the response calculated 
from each band separately. Incident photons can create 
anisotropic charge fluctuations according to the direction 
of the polarization light vector, and those charge fluctu- 
ations relax by emitting a scattered photon and redis- 
tributing charge density via intrinsic electron scattering 
mechanisms such as electron-impurity, electron-phonon, 
electron-electron interactions, or via breaking of Cooper 
pairs. The anisotropy of the charge fluctuations created 
with light can be controlled by aligning incident and scat- 
tered photon polarization vectors, transforming accord- 
ing to the elements of the irreducible point group of the 
crystal. For a material with D 4h tetragonal symmetry, 
and a single Fermi sheet, the B\ g and B 2g Raman re- 
sponses are not coupled to the long-range Coulomb inter- 
action. As a consequence, the Raman charge densities for 
these channels do not couple to the pure charge density 
channel, and the terms given in Eq. (8) vanish. However, 
A\ g fluctuations need not vanish over the unit cell, and 
therefore they can couple to isotropic charge density, giv- 
ing the finite backflow represented by the second term in 
Eq. (7) . We will see that this can also occur for the B\ g 
case for Raman vertices which are allowed when one has 
multiple Fermi surfaces as formed in the Fe-pnictides. 

From the expression for the Raman response we see 
that in the case of an isotropic gap, Ak = A, there 
should always be a peak at 2A. In an unconventional 
superconductor, depending on the polarization, this ab- 
sorption peak is replaced by a peak or other structure 
at 2A , twice the maximum of the gap over the Fermi 
surface. The peak will be rounded by scattering, but 
still provides a measure of the magnitude of the gaps in 
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FIG. 2: The anisotropic energy gaps around the 4 Fermi sur- 
face sheets as a function of angle 6 shown in Fig. 1. For the 
f3 sheets, the solid line is for rp = 1.4 and the dotted line 
r-0 = 0.6 (see Eq. 14). 



the system and may be compared to those determined 
from other experiments, e.g. ARPES and tunnelling. 
The B\g and Bi g vertices in Eq. (3) have zeros in It- 
space and therefore weight the part of the Brillouin zone 
away from these zeros. This is observed, e.g. in the d- 
wave cuprates, where the sharp 2A peak occurs in the 
Bi g channel only, while a less pronounced feature corre- 
sponding to a change in slope occurs at the same energy 
in the Bi g channel which weights the nodal regions most 
strongly. Furthermore, the existence of nodes in a gap 
creates low energy quasiparticles that cause a nonzero 
response for all frequencies. This is in sharp contrast 
to fully gapped superconductors whose response show a 
sharp gap edge with no low energy quasiparticles. 

In what follows we consider separate cases of increasing 
complexity in order to display what features for Raman 
scattering in a multi-band system with different gap sym- 
metries might be expected generically. 



B. Results for some special cases 

1. Single Fermi sheet 

In order to understand the type of behavior found for 
the multiband models of the Fe-pnictides, it is useful to 
begin with some special cases. Early on, motivated by the 
proximity of the Fe-pnictides to a (tt, 0) spin density wave 
phase and the multiple Fermi surface structure found in 
LDA (local density approximation) calculations, a sign 
reversed s-wave gap was proposed. 35 It was suggested 
that spin-fluctuation scattering of electron pairs between 
the a and (3 fermi sheets could lead to pairing with an 
isotropic s-wave gap Ao that changed sign between the a 
and (3 Fermi surfaces. In this case, the Raman response 
would consist of a peak onsetting at 2Aq or several peaks 




FIG. 3: The Raman response of the state with energy gap 
Eq. 10 for r = 1.4. Black/red lines denote B\ g /B2 S , and 
green/blue is the unscreened/screened Ai g , respectively. Note 
that our symmetry classifications are according to the geom- 
etry shown in Fig. 1. For lattice classifications, B\ g and B-^g 
should be interchanged. 



if there were gaps of different magnitudes on the various 
Fermi surfaces. In the absence of impurity scattering and 
inelastic lifetime effects, this peak would vary as (lu — 
2Ao) -1 / 2 as u> approaches 2A Q from above. 

Alternatively, RPA spin-fluctuation calculations and 
functional renormalization group studies find anisotropic 
s-wave gaps which can even have nodes on the (3 Fermi 
surfaces as well as nearby d-wave gaps with nodes on the 
a Fermi surfaces. Here, in order to examine the Raman 
signatures of such states, we consider the simple param- 
eterization 
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This gap is plotted as a function of angle around a circu- 
lar Fermi surface. Note, however, that for a single Fermi 
surface centered at F, this state would not have four-fold 
symmetry. Instead, the Fermi surface parameterized by 
9 is intended to represent a model for the (3i sheet of 
the pnictides, and the Fermi surface angle 9 is measured 
around (tt, 0) rather than T. When combined with the 
/?2 sheet at (0,7r), the full s-wave symmetry of the state 
is restored. The gap in Eq. 10 is normalized such that 
Ao is the maximum over the Fermi surface, and plotted 
for the values of r — rp shown in Fig. 2. For r > 1 
the state has nodes on the Fermi surface, and for small 
values of r — 1 these nodes move towards and 7T. For 
r < 1, there are no nodes but for any nonzero r one has 
an anisotropic gap. 

In Figs. 3 and 4 we exhibit the Bi g , Bi gi unscreened 
Ai g , and screened A± g Raman responses for the gaps 
given by Eq. (10), which are shown in the leftmost 
panel in Fig 2. In Fig. 3, where r = 1.4, there are 
gap nodes on the j3 sheets and one finds the expected 





FIG. 4: The Raman response of the state with energy gap Eq. 
10 for r = 0.6. Black/red lines denote Bi g /B2 g , respectively, 
green/blue is unscreened/screened A\ g . 



low frequency power law behavior in which both response 
functions vary as u>, following the low energy behavior of 
the density of states, since the nodes of the energy gap 
do not align with the nodes of the Raman vertices. For 
the gap edge to = 2A max = 2Aq, the B\ g and unscreened 
Ai g spectra have a log | u> — 2A | singularity. For the 
Big response, there is a change of slope at 2Ao since the 
nodes of the vertex align with the gap maxima. A sec- 
ondary log singularity appears at u> = 2A m ;„, probed 
by the B lg and A lg vertices, but not B 2g - Screening, as 
pointed out before, removes all log singularities from the 
unscreened A\ g response, leaving only a cusp-like behav- 
ior at uj = 2A m i n , 2A max . 

For a nodeless anisotropic gap (r = 0.6) shown as the 
dashed line in the (3\ panel of Fig. 2, the Raman spectra 
are shown in Fig. 4. In this case, there is a gap in the low 
frequency spectra. Here, for r < 1, A m ;„ occurs where 
the magnitude of the gap has a local minimum rather 
than a local maximum as it does for r > 1. This leads 
to a step discontinuity at u>/Aq = 0.5 for the B\ g and 
Ai g spectra rather than the log singularity seen in Fig. 3 
for r = 1.6. The nodes of the Raman vertex for the B 2g 
response are aligned with A m i n , leading to a linear onset 
at uj = 2A min rather than a step onset. 



2. Two Fermi surface sheets 

As indicated by Eq. 4, the A\ g Raman response for 
the case of a multi-sheet Fermi surface is not in general 
additive due to the Coulomb interactions and associated 
charge backflow. While no such backflow appears in the 
Bi g or Big channels for the present case, we will see in 
Sec. Ill that for Raman vertices which are allowed for the 
Fe-pnictides, there can be Coulomb backflow contribu- 
tions to the Big channel. To illustrate this we calculate 
the Aig Raman response for two Fermi surface sheets, 



FIG. 5: The screened A\ g Raman response for 2 bands with 
Ai, 2 (6») proportional to (1 + r cos(46»))/(l + r) with r = 0.2, 
and the maximum gap on Fermi surface 1 equal to Ao and 
on Fermi surface 2 equal to Ao/2. ( black = screened A\ g for 
band 1, red = screened A\ g for band 2, green = mixing term, 
blue = total response.) 



where the energy gaps Ai j2 (0) are both proportional to 
(1 + r cos(40))/(l + r), with r = 0.2. For simplicity, we 
take the gap maximum on one Fermi surface to be Ao and 
on the other Ao/2. We also take equal Raman vertices 
and density of states for each band. 

Expanding Eqs. (4-9) for the case of 2 bands, the over- 
all Raman response can be written as 
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Thus for Ai gi the screened Raman response can be con- 
sidered as a sum of the screened response for each band, 
plus a mixing term A%. Here one can see that if the en- 
ergy gaps and Raman vertices are the same for each band 
the mixing term vanishes, while for all other cases it is 
finite, reflecting the contribution from charge backflow. 
In this case, the light scattering induces interband charge 
transfer in order to recover particle number conservation 
when the charge fluctuations differ on the two bands. 

As illustrated in Fig. 5, the screened Ai g Raman spec- 
trum consists of contributions from each Fermi surface 
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FIG. 6: The Bi g ,B2 S and screened Ai g Raman spectra for 
a four-sheeted Fermi surface, with rp = 1.4. The black line 
corresponds to each contribution from the /3i t 2 bands, the 
red line to the ai band, the green line a<2, and the blue line 
is the total response. Note that for the Ai g case (bottom 
panel) there is an interference term in addition to the sum of 
the contributions from each band. The energy gaps used are 
defined in the text. 



with a gap scale differing by a factor of 2 and an interfer- 
ence contribution coming from the charge backflow. All 
singularities associated with the values of the gaps at the 
stationary points that would appear in the unscreened 
Ai g channel are removed and replaced by cusp-like be- 
havior as in the single band case. However the structure 
of x" c ( w ) i s changed by the interference. 

Interference terms also occur if the gaps are identical 
on each sheet but the Raman vertices differ. A more 
detailed examination of the role of charge backflow is 
presented in Ref. 27. 



III. THE FOUR FERMI SURFACE RAMAN 
SPECTRA 

We would now like to extend this discussion of the 
features in the Raman intensity to the 4-Fermi surface 
model of the Fe-pnictides. Here we first consider the two 
sets of A„(0) gap variations for the four Fermi surface 
sheets v = ai,a 2 ,f3\ and f3 2 shown in Fig. 4. Specifically, 
the energy gaps are taken to be (Fig. 2) 



A /3i,ft! = A 



A Ql = -0.8A, 



1 T rp cos(20) 

1+773 

1 + cos(4#) 



A, 



-0.4A, 



1 - cos(46>) 



(14) 



As previously discussed, these gap choices are motivated 
by the anisotropic gaps found in RPA spin fluctuation 
and functional renormalization group calculations. The 
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FIG. 7: Same as Fig. 6 but for rp = 0.6. 



amplitudes of the a\ and a 2 gaps have been chosen to 
avoid an accidental overlap of singularities between the a 
and (3 gap extrema. The first of these, shown by the solid 
curves in Fig. 2, corresponds to an A\ g gap with nodes 
on the j3- Fermi surfaces (rp = 1.4), and the second one is 
nodeless corresponding to the dashed curves (rp = 0.6). 
The Bi g ,B 2g and A\ g spectra for these two cases are 
shown in Figs. 6 and 7, respectively. The contribution 
from the individual Fermi surfaces are also indicated. 

For the B lg and B 2g spectra the structure seen in the 
total response are just the sum of the spectra from the in- 
dividual Fermi surface sheets with the appropriate gaps 
shown in Fig. 2. For example, in Fig. 6, the B\ g re- 
sponse for the gap with nodes consists of a sum over 
different contributions coming from the two hole Fermi 
surfaces ol\ and a 2 and the sum of the spectra from the 
Pi and (3 2 electron Fermi surface sheets which are iden- 
tical since Raman only probes | A | and is not sensitive 
to the phase in the absence of impurities. Just as the for 
the previous discussion of the single Fermi surface case, 
one can easily identify the characteristic features com- 
ing from each Fermi surface shown in Fig. 6. The B\ g 
Raman response for a gap with nodes on the (3 sheets 
(red curve) exhibits log singularities at u / Aq equal to 
2Ap, max /A ~ 1.5 and 2 | A^ \ min /A ~ 0.25. The ai 
sheet (green) contributes an additional log singularity at 
2 | A ai Imaa; /A ~ 1.6. Since the B\ g vertex has a node 
at the minimum value of the gap on the ol\ sheet the 
discontinuity at ioj Aq = 0.8 is eliminated leaving only 
a linear onset. The contribution of the a2-sheet has a 
similar linear onset at u> = and a change in the slope 
at 2 I A 



Q2 \max 



/A ~ 0.8 due to the fact that the B\ g 
Raman vertex has a node at the minimum and maximum 
values of the gap on the a 2 Fermi surface. 

A similar discussion can be given for the B 2g case 
shown in the center panel of Fig. 6, where the total re- 
sponse is the sum of the contribution from the individual 
Fermi surface sheets with the B 2g Raman vertex. As dis- 
cussed in the previous section, the screened A± g response 
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FIG. 8: Same as Fig. 6 but for different Raman vertices, given 
by Eq. (15). 



shown in the lower panel of Fig. 6 contains the screened 
contributions from each band, plus an interference term 
with interference contributions coming from each pair of 
bands. All singularities are removed by the backflow, 
leaving the cusp-like behavior which can be tied to each 
of the gap extrema, as done for the other channels. 

As mentioned in Sec. II A, different forms for the Ra- 
man vertices are in principle allowable other than the 
ones given in Eq. (3), giving the Raman spectra shown 
in Figs. 6 and 7. Here we consider how different forms 
for the Raman vertices affect the general structure of the 
channel dependent Raman spectra. As an example, we 
consider the following form for the Raman vertices: 
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= -l, 




= cos(26>) 




= ±1, 




= sin(20) 




= sin(0). 



(15) 



These vertices are all allowed by symmetry and have 
less anisotropy than the ones considered in Eq. 3 and 
thus highlight different gap structures. Moreover, due to 
the angle-independent form for the B\ g vertices on the (3 
sheets, screening must be included and the mixing terms 
are non-zero to give finite contributions to the spectra. 

The resulting channel-dependent spectra are shown in 
Fig. 8 for the Fermi sheet gaps given in Eq. 14 with 
rp = 1.4. For the B\ g case, Raman scattering from the 
individual (3 sheets now is canceled by backflow due to 
the angular independent Raman vertices on those sheets, 
and two identical contributions arise from the a sheets, 
as in Fig. 6. However, an additional mixing term aris- 
ing from scattering interferences involving each separate 
(3 sheet, gives a large contribution to the spectra, with 
peaks occuring at the energy scales given by the energy 



FIG. 9: Polarization dependence Raman response for the s± 
state, having energy gaps given by Eq. (16) and Raman ver- 
tices given by Eq. (15). Color scales, denoting contributions 
from each band and the total response, are the same as in 
Figs. 6-8. 



gap extrema on each sheet. For the Bi g spectra, the si- 
nusoidal variation sin(0) of the vertices on the f3 sheets 
now allow for the highlighting of the gap maximum, giv- 
ing a peak frequency at twice the gap maximum for the 
j3 sheets, in contrast to the spectra shown in Fig. 6. A 
dramatic change of the spectra is observed for the A\ g 
channel. Due to backflow, the Raman response from each 
separate band vanishes identically for the constant ver- 
tices in Eq. (15), but due to the change in sign, the entire 
Raman spectra arises solely from the mixing terms. The 
positions of the cusp- like feature at lu — 2A are the same 
as that shown in Fig. 6 for a different form for the Raman 
vertices, however the overall spectra lineshape is qualita- 
tively different. The overall structure of the lineshapes 
thus indicate that a proper account of the anisotropy of 
the Raman vertices may be needed in order to obtain a 
qualitative comparison with the experimental observed 
Raman spectra. This is a topic for further study. 

Finally we consider a simple sign-changing s± state on 
the four sheets: 



= Ao, 



= -Ao/V2, 



A, 



A /4. 



(16) 



Using the vertices defined by Eq. 15, the resulting Ra- 
man spectra is shown in Fig. 9. Here, the square root di- 
vergence at twice the gap value is ubiquitous, displaying 
both in Bi g and in Bi g The only polarization difference is 
that the contributions from the (3 bands for the B\ g chan- 
nel are screened, in contrast to the B\ g channel. Apart 
from this difference, the spectra are qualitatively similar, 
due to the mixing terms in B\ g which restore the singu- 
larity at u> = 1.5A . For the Ai g channel, the divergences 
are screened and the response maintains thresholds and 
peaks at twice the gap energy for each band. 
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IV. CONCLUSIONS 

Here we have studied the Raman scattering response 
for some simple models of the Fe-pnictidc superconduc- 
tors. Specifically, we considered two different anisotropic 
A\ g gaps, one with nodes and one without nodes, on 
four circular Fermi surfaces. Besides the well-known low 
frequency differences in the spectra for nodal and non- 
nodal gaps, we found a rich set of high frequency struc- 
tures arising from stationary points of (0) on the var- 
ious Fermi surfaces. Measurements of different polariza- 
tions may allow one to associate particular gap structures 
with individual Fermi surface sheets. If the gap has a 
large anisotropy as suggested by some calculations, there 
will be a rich Raman spectrum for different symmetry 
channels. However, if one has a relatively isotropic gap, 
such as the proposed sign-switched s— wave, the spectrum 
should be simpler and sharper since for an isotropic gap 
the response has a square root singularity at 2Aq rather 
than the weaker log singularity found for an anisotropic 
gap. 



We have also discussed some of the unusual aspects of 
the Raman spectra to be anticipated in the Fe-pnictides 
due to their multisheeted Fermi surface electronic struc- 
ture. In particular, Coulomb backflow mixing may affect 
the Big spectrum as well as A\ g . There are other un- 
usual aspects, such as excitonic modes associated with 
short range interactions 38 or competing order parame- 
ter channels, which we have not explored here, but may 
make the spectra in these materials even richer. Work 
along these lines is in progress. 
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